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Zusammenfassung 



Axionen entstehen bei der Losung des „starken CP-Problems" mit Hilfe der Peccei- 
Quinn-Symmetrie. Der astrophysikalisch und kosmologisch prinzipiell erlaubte Bereich 
von Axionenmassen beschrankt sich auf zwei Regionen: zum einen „leichte" Axionen mit 
Massen zwischen 10 -5 und 10 -3 eV, die kalte dunkle Materie (CDM) des Universums 
darstellen, zum anderen „schwere" Axionen aus dem sogenannten „hadronischen Axio- 
nenfenster" mit Massen um 10 eV. Ihre Existenz wiirde mit £l a ~ 0.1-0.3 zur heifien 
dunklen Materie (HDM) beitragen. 

Die genaue Breite des hadronischen Axionenfensters hangt stark vom zugrundelie- 
genden Axionenmodell ab. Deswegen werden in dieser Arbeit zunachst die bestehenden 
und fur hadronische Axionen relevanten astrophysikalischen Argumente neu uberpriift. 
Das Ergebnis ist, dass Axionen mit einer Masse zwischen 10 und 20 eV nicht aus- 
geschlossen werden konnen, falls sie keine Baum-Niveau-Kopplungen an gewohnliche 
Quarks und Leptonen besitzen und wenn gleichzeitig ihre Kopplung an Photonen ext- 
rem stark unterdriickt ist. 

Hadronische Axionen wechselwirken in erster Linie mit Protonen und Neutronen. 
Als Konsequenz konnen sie in nukleonenreicher Materie durch den Bremsstrahlungspro- 
zess N + N —* N + N + a erzeugt werden. Daher - um das hadronische Axionenfenster 
genauer zu untersuchen - betrachten wir Akkretionsscheiben um Schwarze Locher 
und isolierte Neutronensterne. Erstere konnen wahrend des Kollaps eines Neutronen- 
doppelsterns entstehen und liefern ein vielversprechendes Modell fur kurze Gamma- 
strahlungsausbruche. Gemafi diesem entstehen in der Akkretionsscheibe Neutrinos und 
Antineutrinos, von denen sich ein Teil in Elektron-Positron-Paare und weiter in hoch- 
energetische Gammastrahlung umwandelt. Ein zusatzlicher Energieverlust durch die 
Emission hadronischer Axionen konnte die Luminositat der Gammastrahlungsausbruche 
erheblich verringern. Wir stellen fest, dass Axionen mit einer Masse von etwa 11 bis 
66 eV einen wichtigen Einfiuss auf Gammastrahlungsausbruche haben, wobei diese 
quantitativen Aussagen jedoch stark vom zugrundeliegenden Modell der Akkretions- 
scheibe und dem des Schwarzen Lochs abhangen. 

Signifikantere Grenzen fur die Masse hadronischer Axionen ergeben sich aus der 
Betrachtung isolierter Neutronensterne. Mit einem einfachen Modell fur das Abkiih- 
lungsverhalten eines Neutronensterns zeigen wir, dass die gemessenen Oberflachentem- 
peraturen der Pulsare PSR 1055-52, PSR 0630+178 (Gemina) und PSR 0656+14 zu 
hoch sind, um mit der Existenz hadronischer Axionen im Einklang zu stehen. Diese Er- 
kenntnis wird durch numerische Berechnungen der Neutronensternkiihlung von Umeda 
et al. H| bestatigt, deren Ergebnisse wir einer erweiterten Interpretation im Hinblick 
auf hadronische Axionen unterziehen. 
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Chapter 1 
Introduction 



The axion was proposed in 1977 to solve the CP problem of QCD ("strong CP 
problem"). Since then, more than 500 research papers about axions and related 
topics have been published 0. Evidently, axions have been a stimulating subject! 
The main reason for this steady interest is that axions not only provide an elegant 
solution to the strong CP problem, but are also attractive from a cosmological 
perspective. More specifically, there are two regions of axion masses where these 
particles could constitute some or all of the cosmic dark matter: Axions with 
a mass in the range m a pa 10 _5 -10~ 3 eV are an ideal candidate for cold dark 
matter (CDM), and two current experiments actually search for these very light 
particles |2], ||. Moreover, it was recently claimed that "heavy" axions in the so- 
called "hadronic axion window" around 10 eV provide a hot dark matter (HDM) 
component with Q a pa 0.1-0.3, exactly the amount needed in mixed dark matter 
scenarios ||. 

While there is no doubt that axions in the CDM window are astrophysically 
and cosmologically allowed f|, the situation is less clear for the HDM range 
because the existence and exact width of the hadronic axion gap strongly depends 
on details of the axion model. Therefore, firstly we have re-examined the hadronic 
axion window and have found that axions without tree-level couplings to ordinary 
quarks and leptons ("KSVZ-type axions") and masses between 10 and 20 eV are 
not excluded by previous astrophysical arguments if their coupling to photons g ail 
is less than approximately 3 x 10~ n GeV -1 . This value is nearly one hundred 
times smaller than in generic axion models, implying a strong suppression of the 
photon coupling for these models to be viable in the hadronic mass range. It turns 
out that one can construct models where this suppression arises as a accidental 
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cancellation between two independent contributions to the axion-photon coupling 
|J. Therefore the hadronic axion window is indeed open. Henceforth, the term 
"hadronic axions" shall refer to models with severely suppressed photon couplings, 
i.e. they are a subclass of "KSVZ-type axions," a term which applies to models 
without tree-level coupling to ordinary quarks and leptons. 

Hadronic axions are very weakly interacting particles, and one might think 
that they can not have any significant influence on astrophysical objects. How- 
ever, as a consequence of the axion-pion mixing, their effective nucleon coupling 
is comparable to those of generic axion models. This means that those astro- 
physical objects which mainly consist of nuclear matter are suitable candidates 
to explore the hadronic axion window. Therefore, we consider isolated neutron 
stars and black hole accretion discs (BHADs) as axion laboratories — both types 
of objects have not been studied in this context. 

Accretion discs around black holes are possible outcomes of neutron star merg- 
ers and have recently been discussed in connection with a very intriguing class of 
mysterious astrophysical phenomena: the gamma-ray bursts, short and intense 
photon eruptions with energies of typically 100 keV-1 MeV. According to the 
BHAD model, a gamma-ray burst develops when neutrinos, which are produced 
inside the hot accretion torus, are converted into a photon burst. The crucial 
point is that the initial neutrino luminosity must be extraordinarily large to ac- 
count for the observed energies in photons in the range 10 50 -10 53 erg. However, 
numerical neutron star merger simulations |7j yield total neutrino luminosities 
of about 10 53 erg s _1 , enough to explain short and weak gamma-ray bursts with 
energies in the ballpark of 10 50 -10 51 erg. This situation could change if axions 
existed in the hadronic axion window. Produced via the bremsstrahlung process 
N + N^N + N + a, they represent a novel, important energy-loss mechanism 
in accretion discs. In this case, the neutrino luminosity would be significantly 
reduced, dimming the observed gamma-ray burst signal. Put another way, the 
requirement that the axion luminosity must not exceed the neutrino emission puts 
an upper axion mass limit, which we find in the range 10-100 eV, i.e. hadronic 
axions could have an impact on the evolution of gamma-ray bursts. However, 
this bound depends on unknown details of the BHAD model. Nevertheless, al- 
though precise and reliable statements about the influence of hadronic axions on 
gamma-ray bursts require a deeper understanding of these fascinating celestial 
objects, it is quite certain that hadronic axions if they were found have to be 
included in present BHAD models. 
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Finally, we discuss the cooling of neutron stars by means of hadronic axion 
emission. We set up a simple cooling model and show that recently measured 
surface temperatures of the isolated neutron stars PSR 1055-52, PSR 0630+178 
(Gemina), and PSR 0656+14 are too high to be in accordance with hadronic 
axion emission because these particles significantly accelerate the energy loss of 
the star. Furthermore, we review a recent calculation of Umeda et al. || who 
studied neutron star cooling including axions by means of numerical simulations. 
However, they did not discuss the implications for hadronic axions. It turns out 
that an interpretation of their results with regard to the hadronic axion window 
confirms the outcome of our simple model, namely that hadronic axion emission 
is in conflict with the observed surface temperatures of neutron stars. 
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Chapter 2 

Summary of Axion Physics 



By way of introduction, we briefly summarize the most important aspects of axion 
physics. First, we outline how axions solve the strong CP problem, and how they 
arise as a simple extension of the Standard Model. Then, the KSVZ axion is 
discussed in more detail. Finally, the couplings of axions to ordinary matter, i.e. 
to photons, electrons and nucleons are discussed. 



2.1 Motivation for Axions: The Strong 
CP Problem 

In weak interactions one observes CP violation, whereas the strong interactions 
are known to respect P and CP invariance to very high accuracy. This CP 
invariance of the strong sector is a problem because QCD predicts CP violation 
as a consequence of two different mechanisms. First, there is the topologically 
nontrivial ground state of QCD, the so-called "0- vacuum". Second, the axial 
transformation, necessary to diagonalize the complex quark mass matrix, also 
leads to CP violation. Both effects together cause an additional nonperturbative 
term in the QCD Lagrangian, 

n 2 

*-SM — ^SM 327T 2 ^' 

where £g^f is the perturbative part of the Standard Model Lagrangian, g s the 
strong coupling constant, G% v the color field strength tensor, G a = ^e M „ pCT G af " T 
its dual, and = 6 + Arg detM. The e tensor that occurs in the definition of 
G implies that is not invariant under parity, and thus odd under CP. Note 
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that the effective CP violating parameter G is the sum of a QCD contribution — 
the vacuum angle — and an electroweak part — Arg det M — related to the phase 
structure of the quark mass matrix M. The experimental bound on the neutron 
electric dipole moment requires that 

|9 + Arg detM|<l(r 9 . (2.2) 

Now, it is difficult to understand why a compound quantity like 9, which is a sum 
of two very different contributions, should be so small. This fine-tuning problem 
is known as the strong CP problem. 

The strong CP problem can be elegantly explained by introducing an addi- 
tional global, chiral symmetry with an associated current and charge: the Peccei- 
Quinn-symmetry U(1)pq and the Peccei-Quinn charge ||. Its existence requires 
certain enhancements of the Higgs sector of the Standard Model so that its en- 
larged symmetry group becomes SU(3)c x SU{2) L x U(l)y x {7(1)pq. However, 
it should be stressed that C/(1)pq is a global symmetry and not a gauge symme- 
try. In typical axion models, PQ symmetry is achieved by introducing additional 
Higgs fields $ with degenerate vacua; at least one is necessary. However, in the 
real world the chiral symmetry C/(1)pq is not observed, hence it must be sponta- 
neously broken at an energy scale f a ; the associated Nambu-Goldstone boson is 
the axion. Under a PQ transformation the axion field a(x) shifts as 

a(x) — > a(x) + a f a , (2.3) 

where a is the parameter associated with the U(1)pq transformation. Now we 
come to the essential feature of the PQ mechanism: the PQ symmetry U(1)pq is 
only exact at the classical level. At the quantum level, however, it is broken by 
the Adler-Bell-Jackiw anomaly because the PQ symmetry is chiral. As a result, 



an additional term appears in the effective Lagrangian of equation (pTT 



= ^sm + e GT G% + £ ^32^2 °a % » (2-4) 

where £ is a model dependent parameter. The £-term is very similar to the CP 
violating B-term and provides an effective potential V e s(a) for the axion field a(x). 
This potential implies that the degeneracy of the vacuum is removed. Peccei and 
Quinn showed that the vacuum expectation value of the axion field, i.e. the value 
of a(x) where V e s(a) has its minimum, is given by 

(a(x)) = -|e. (2.5) 
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Now one introduces a physical axion field a p h y , which is defined as excitations 
around the vacuum expectation value (2.5), i.e. 



a{x) phy = a(x) - (a(x)) , (2.6) 



and substitutes it into expression ( |2.4| ). In doing this, the QGG term is elimi- 
nated, and therefore the strong CP problem is solved. Furthermore, one is left 
with an axion-gluon interaction of the form a^yGG — obviously an immediate 
outcome of every axion model. As a result of this generic axion-gluon coupling, 
axions pick up a small mass of order 

m a « — — 

la 

Ar AO 6 GeV\ , „, 

« 6 eV ( — J , (2.7) 

where f n ~ 93 MeV is the pion decay constant and = 135 MeV the pion mass. 
As the axion has a mass, it is actually a "pseudo Nambu-Goldstone boson." 

In summary, we arrive at the following Lagrangian, where from now on a 
stands for the physical axion field a p h y , 

<eff - ' 



£sm = £sm + -d tl ad t "a 



d^a \ a g 



2 



+ M±'^) + & GirG - (2 - 8 » 

The second term is the kinetic energy of the axion field, the third represents 
possible interactions of the axion with fermions ip and has to be in agreement 
with the classical invariance under {7(1)pq, i.e. under the transformation ( |2.3j ). 
Hence we conclude that the axion field a may enter the interaction Lagrangian 
£i n t only through derivative terms d^a. 



2.2 Axion Models 

2.2.1 Standard Axion and Invisible Axions 

Let us now look at different axion models which are based on these general ideas. 
First, we consider the question how the axial symmetry {7(1)pq, which solves 
the strong CP problem, can be reconciled with the Standard Model. In the 
Standard Model one introduces a complex scalar field, the well-known Higgs field 
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0. It is an SU(2)l doublet, i.e. its weak isospin is Iw = \- Furthermore, its 
weak hypercharge is Y = 1. As a consequence of the Higgs mechanism, all 
boson and fermion masses are generated by this Higgs field 0. The four phases 
of the isodoublet provide four additional degrees of freedom. Three of them 
provide the longitudinal degrees of the W and Z°, and the remaining phase 
ends up as the Higgs boson. So, in the Standard Model there is no room for 
axions. Another way to see this is to recall the transformation properties of the 
Standard-Model Lagrangian under axial transformations. In particular, consider 
the Yukawa-coupling terms 

£ Yu k = ~ G d qL<frd R - G u q L C u R + h.c. , (2.9) 

where q = (u, d) and C = 2T 2 0*. One can ensure that the first term is invariant 
under the special axial {7(1)pq = U(1)r — U(1)l transformation [] q L — > e ta ^ 2 q^, 
dji —> e~ ta l 2 (Ir, ur — > e~ lOL l 2 ur if one requires the Higgs field to transform as 

0^e ia 0. (2.10) 

However, the second term in ( |2.9| ) is only invariant if C transforms in the same 
way as 0, but this is not the case. 

A simple solution is to introduce a second Higgs doublet 2 in place of <j> c , 
which is therefore independent of <p\ = <ft. Then, both Higgs fields <p\ and 02 
transform as in ( p,10| ), and PQ invariance of ( |2.9|) is attained. Finally, when 
the SU(2)l x U(1)y x U(1)pq symmetry breaks spontaneously, the massless ax- 
ion appears together with three other Nambu-Goldstone bosons, which will be 
"eaten" by the massive and Z° gauge bosons. This is the original standard- 
axion model in which the breaking scale of the PQ symmetry is equal to the 
electroweak one, i.e. f a ~ 250 GeV. The associated axion with a mass of about 
100 keV, however, was quickly ruled out by experiments: if weak-scale axions 
were to exist, one would expect decays in the meson system like 

K± -> tt ± + a . 

1 Consider a general axial transformation U(1)a and its related axial vector current j'i — 
?p-f^-f 5 ip. With ipL = 5(1 ^7 5 )V' an( l — 5(1 +7 5 )V' one can write the axial current in terms 
of left- and right-handed currents j£ and j^, 

3a = (ipR + ^l)i^1 5 {^r + i/>l) = tpR^ipR - TpLl^lpL = 3r - 3l> 

where we have used 7 5- 0fl = ipR, J 5 ipL = — and tpLj^ipR. = ipR.J^ipL = 0. Therefore we may 
write U(1)a = U(1)r-U(1)l. 
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The absence of such decays implies that the standard axion can not exist. Even 
worse, its mass must be less than about 10 keV, or equivalently the symmetry- 
breaking scale must be larger than 1000 GeV. 

In order to maintain the PQ solution of the strong CP problem, the "invisible 
axion" was invented. The special feature of these models is that the axion resides 
in the phase of a new complex SU(2)l x U(1)y singlet scalar field $, which has 
a nontrivial PQ charge Q$>. As a consequence, the field $ does not participate in 
weak interactions, and hence the PQ symmetry breaking is decoupled from the 
electroweak one. Therefore, the PQ scale f a is an arbitrary parameter in invisible 
axion models, implying also that the axion's couplings are not fixed. Perhaps the 
most common axion models are those of Dine, Fischler, Srednicki, and Zhitnitskii 
(DFSZ) and Kim, Shifman, Vainshtein, and Zakharov (KSVZ). The KSVZ axion 
is an example of a wider class of axion models, the so-called "KSVZ-type axions." 
The main difference between the DSVZ and KSVZ-type models is that the latter 
have no tree-level couplings to ordinary quarks and leptons. This is achieved by 
setting the PQ charge of all ordinary fermions to zero. Instead, one introduces 
new, heavy, colored fermions with a nonvanishing PQ charge Xf. 

2.2.2 The KSVZ Model 

In order to gain deeper insight into axion physics, one specific example shall be 
discussed in some detail, i.e. the KSVZ model. Since the KSVZ axion has no 
tree-level coupling to ordinary quarks, the question of whether such an axion is 
able to solve the strong CP problem arises: axions need an anomalous triangle 
coupling to two gluons, i.e. a U(1)pq-SU(3)c-SU(3)c anomaly, which is a fun- 
damental feature of every axion model. Therefore, within the KSVZ model, one 
has to introduce at least one new heavy Quark Q with nonvanishing PQ charge 
and nontrivial transformation properties under SU(3)c, e.g. a SU(3)c triplet. 
As a result, the essential aGG coupling can be realized at lowest order via an 
anomalous heavy quark triangle loop, which is illustrated in Fig. |2.1| . 

Let us now consider the spontaneous breaking of the C/(1)pq symmetry in 
the KSVZ model. To this end, we introduce a new field $ with the associated 
Lagrangian 

= (<9 / ,$) t (<^$) + /p ; Q $ t $-A($ t $) 2 , (2.11) 
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Figure 2.1: Interaction of axions with gluons via a heavy quark loop. 

where A, /|q > 0. Below we set A = 1. is invariant under 

$ -> e to $, (2.12) 

i.e. possesses a C/(l) global symmetry, the PQ symmetry £/(1) P q. The poten- 
tial V($) then has a circle of minima at ($) = {fe^/ \f2)e tLp with an arbitrary 
phase v?. Now we spontaneously break the U(1)j>q symmetry by expanding $ 
around one of these ground states ($) with minimum energy. If we introduce 
new fields a and p corresponding to angular and radial excitations around ($), 
we may write 

$ = J=(/p Q + p)expf-^V (2.13) 



\/2 V/pq 
The massless mode a is linked with the axion field and translates under {7(1)pq 
in ( p.!2j ) as a — > a + a/pQ. Now we substitute (|2.13| ) into (|2.11|) , i.e. we express 



the Lagrangian in terms of a and p. The result is 

+7^P(^«)(^ a ) + 7T7^P 2 (^«)(^«) + const. (2.14) 
/PQ 2 /pq 

The third term has the form of a mass term (— |m^p 2 ) for the p field. Thus, 
the p mass is m p = v^/pq- Because /pq is a high-energy scale, p is a very 
heavy particle and can be neglected as far as the low-energy regime is concerned. 
Therefore, we may keep only the axionic degree of freedom, i.e. we drop all terms 
involving p. 

Finally, we have to implement the quark Q into our model. To this end, we 
add a usual Dirac term for the fermionic Q field and introduce Yukawa couplings 
between Q and $. We then write 

C Q = % - (Qi^Q + h.c.) - h(Q L Q R $ + &Q R Q L ) (2.15) 
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with a Yukawa coupling constant h > 0. Remember, the {7(1) P q transformation 
of the $ field is already fixed through ( 2.12 ). Hence, the transformation property 
of the quark field Q is determined by the demand for C/(1)pq invariance of the 
Yukawa interactions ( 2.15 ). This invariance can be achieved by 



Qi 



Ja/2 



Ql, 



Q 



R 



-ia/2 



Qr, 



(2.16) 



which is equivalent to the axial transformation Qr + Ql = Q — ► e~ ir ' a ^ 2 Q. Note 
that both Q and Q transform in the same way, i.e. it is also Q — > Qe" 17 '"' 2 . If 
we insert ( |2.13| ) into ( |2.15| ), we obtain for the Yukawa term, i.e. the last term in 

(HID, 



Yuk 



V2 



h (Q L Q/?e M// PQ + e" M// PQ QrQl 



(2.17) 



where we have neglected all terms involving p. One can simplify this expression 
by using the normal field Q = Qr + Ql instead of the chiral quark fields Ql and 
Qr. Recalling that 



QlQr = Q\(1 + 1 5 )Q and Q R Q L = Q~(l ~1 5 )Q, 



(2.18) 



one finds 

^Yuk 



V2 



V2 



hQ 



hQ 



V/PQ 



1+7 5 



+ e" 



V/PQ 



V/pQ 



+ e" 



V/pQ 



+ 7° 



i -r 

2~ 

. e ia /fpQ . 



Q 

e ~ ia /fpQ 



Q 



= Jj^hQ [cos(a// PQ ) + sin(a// PQ )] Q 
= —ilR h Q g^V/pQ q 

With $n§, (p!5|) , and (p!9|) we finally get the KSVZ Lagrangian 

1 * f, '°hQe^ a/f PQQ. 



(2.19) 



£ K svz = | (Q7^Q + h.c.) + J(,9 M a) 



v/2 



(2.20) 



In order to interpret the last term of this expression, we expand it in powers of 
a /fpQ- The zeroth-order term provides a mass for the quark Q. Higher-order 
terms describe interactions between Q and the axion a. One obtains 



^KSvz = -(Q7^2 + h.c.)+-(^ 



fpQ 

V2 



hQQ + Ant 



'.2.2V. 
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and thus a Q mass of mg = hf-pq^/yl. In addition, the interaction Lagrangian 
is given by 

Ant = -^aQ 7 5 Q + £f a'QQ + .... (2.22) 
/PQ z Jpq 

One then can see that to lowest order the interaction between the axion a and the 
exotic quark Q is pseudoscalar. The corresponding Yukawa coupling constant is 
g a Q = tuq/ /pq. There are, however, higher-order terms, which sometimes must 
be taken into account in order to achieve the correct result. Therefore, it is better 
to search for an alternative approach to the a-Q interaction. One may remove 
the axion field in the Yukawa coupling term (last term in fl2.20|) ) by a local chiral 
transformation of the Q field, 

Q e -nV2/ PQ q ; q e -; 7 V2/ PQ q _ ^.23) 

With = e~^ a Y it is straightforward to show that under (|2~23| ) the KSVZ 

Lagrangian (|2.20 ) transforms as 

>Cksvz - ^Q7 M 7 5 Q^a+ l -{d^af - m Q QQ . (2.24) 

Obviously, the transformation (|2.23|) generates a derivative axion interaction from 
the kinetic Q term in £rsvz ( ^.20[ ), namely 

C iBt =-^-Qf^Qd lt a. (2.25) 

In contrast with (|2.22|) , this interaction contains no higher-order terms. The 
existence of two different couplings 

A„t = -^aQ 7 5 Q and C mt = -l-Q^Qd^a (2.26) 

JPQ Z JPQ 

raises the question of which is to use in concrete calculations. It would be naive to 
suppose that both interaction Lagrangians always yield the same result because 
the pseudoscalar one is only the leading term of an expansion in powers of a//pQ. 
Therefore, it is a safe strategy to use the derivative coupling in all calculations. 

To conclude this section we finally consider the effects associated with the 
chiral anomaly, i.e. the axion mass. We go back to the Lagrangian fl2.20|) which 
is invariant under global £/(1)pq transformations ( j2.16 ) and a —>■ a + /pqo:. The 
associated current due to Noether's theorem is simply 

Jpq = /pq^-^27VQ (2-27) 
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0„7'Sn = . (2.28) 



with 

However, the isosinglet axial current |Q7 M 7 5 Q has an Adler-Bell-Jackiw anomaly 
and is therefore not conserved at the quantum level. For the divergence of jp Q 
one finds explicitly [|K| 

32tt 2 



d ^Q = ~4^ G a U G%- (2-29) 



Recall that at the scale /pq the PQ symmetry breaks spontaneously and the 
axion arises as the associated Nambu-Goldstone boson. As a result, the axion 
is a massless particle. However, at energies below Aqcd, the axion develops a 
mass due to QCD effects: the heavy quark Q has a color anomaly. So, the 
axion interacts with gluons (Fig. |2.1|) , and therefore with quark-antiquark states. 
Furthermore, the axion mass should vanish in the limit of vanishing u- or ci-quark 
masses.^ With this information, Bardeen and Tye []TT| constructed an anomaly 
free current jpq out of jpq, 

Jpq = faPa - \q^1 5 Q+ \ ( ^ vrf^u + ™" d-f-fd) , (2.30) 
w 2 2\m u + m d m u + m d J 

where g7 M 7 5 g, q = u,d are the chiral quark currents. In the case that at least 
one quark mass vanishes, the last term of expression ( |2.30| ) is conserved up to the 
chiral anomaly, i.e. 



l( md -uY7 5 u+ mu dYl 5 d 



2 \m u + m d m u + m d 



^G?G% (2-31) 



for m u — > or m d — > 0. Therefore, with fl2.29Q , ( |2.30| ), and ( |2.31| ), it is easy to 
see that the Bardeen- Tye current is anomaly free and thus conserved in the limit 
of vanishing u- and d- quark masses, 

^aJpq if m u o r m d 0. (2.32) 

That means, the conservation law of the Bardeen- Tye current jpq is explicitly 
broken by nonvanishing u- and <i-quark masses. This fact is referred to as the 
partial conservation of the axial current (PCAC). Using PCAC and the standard- 
current algebra approach |I2|, in which the axion mass is related to the "soft" 



2 Effects due to the s-quark are of order w = m u /m s w 0.029 and are neglected for the sake 
of simplicity. 
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divergence of jpq, the axion mass can finally be calculated as 



m a = — —— , (2.33) 

1 ~~ t - Z 



where z = m u /md, f n ~ 93 MeV the pion decay constant and = 135 MeV 
the pion mass. 



2.3 General Axion Couplings 

For the most part, axion physics is determined by the scale f a of PQ symmetry 
breaking. According to ( |2.8|) , the axion's couplings to photons and fermions are 
all proportional to the inverse PQ symmetry breaking scale or, equivalently, to 
the axion mass 

Axion- Couplings ~ — oc m a . 

J a 

In detail, however, these couplings are model dependent, i.e. they depend on the 
implementation of the PQ mechanism. 



2.3.1 Photons 



First we consider the axion's coupling to photons. An apparent reason for this 
interaction is the axion-pion mixing: As a result of the electromagnetic anomaly, 
pions couple to two photons, causing an effective axion-photon coupling. Over 
and above that, there is a second contribution to the axion-photon interaction. 
In general, the axion has Yukawa couplings to fermions which carry PQ charges. 
If these PQ fermions also carry electric charges, they couple the axion to two 
photons by means of a triangle loop (Fig. 2.2). As a consequence of the axion- 
photon coupling, axions decay into two photons with a lifetime 



T, 



£177 



2 3 

a m a 



C 2 

2567r 3 / a 2 a77 




'V\/\/V~V"' 



Vww 



V 



+ 



Figure 2.2: Coupling of the axion to two photons. The left triangle loop is made 
up of fermions carrying PQ and electric charge 
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= 3.53 x 10 24 sec m e " v 5 , (2.34) 

where a = 1/137 is the fine-structure constant. 

Let us now study some characteristic properties of the axion-photon coupling. 
If we demand CP-invariance, the effective Lagrangian for the axion-photon inter- 
action can be written as 

£(177 = ~^9a"/"/F^ v F^ (X Qa-yy & E ■ B , (2.35) 

where F is the electromagnetic field strength tensor, F its dual, a the pseu- 
doscalar axion field, and g ail the axion-photon coupling strength with the di- 
mension (energy) . The latter is given by 

E 2(4 + z + w)' 



a 

Qarff 



N 3(1 + z + w) 



(2.36) 



where a = 1/137 is the fine structure constant, E the electromagnetic anomaly, 
and N the color anomaly of the PQ symmetry. They are given by 

N^Y^Xj, E = 2j2x j Q) Dj , (2.37) 

3 5 

where Xj and Qj are the PQ charges and the electric charges of the PQ fermions, 
respectively, and Dj = 1 for color singlets (charged leptons) and 3 for color triplets 
(quarks). Finally, z and w are the mass ratios of the u- to the <i-quark and the 
u- to the s-quark, respectively. The first term of expression (|2.36|) corresponds 
to the electromagnetic anomaly of the PQ fermions, whereas the second term is 
associated with the axion-pion mixing. 

From ( |2.36| ) it is clear that the light quark mass ratios z and w are of great 



interest, since they determine the axion-photon coupling strength g ail . However, 
even though there has been considerable effort to calculate the numerical values 
for z and w, the results are still controversial and remain under active discus- 



sion. Leutwyler QT3j| used chiral perturbation theory results for the kaon and 
pion masses and extracted from those the relative sizes of m u , m^, and m s . He 
obtained the very stringent constraints 

z = m u /m d ~ 0.553 ±0.043 

w = m u /m s w 0.029 ±0.003. (2.38) 

However, alternative approaches have been made, e.g. using sum rules and numer- 
ical simulations of QCD on a lattice. These calculations lead to different results 



15 



which are not mutually consistent. Therefore, it remains to be seen whether the 
results ( p. 38 ) are reliable. The Particle Data Group gives the conservative range 



z = 0.2-0.7 \TA\. Therefore, we rather use 

z ^0.55 ±0.1 (2.39) 



than the z- values of expression (|2.38|) . With (|2.7|) and ( |2.39|) , we finally obtain 



the axion-photon coupling in terms of the axion mass, 

= 2^(|- L93±a08 

= m>v(^ l93±o 4 (2 - 40) 

where have used m e y = m a /eV and ( |2.38|) . Hereafter we use the abbreviation 



C a77 = ^ - 1.93 . (2.41) 

In the DSVZ model and grand unified theory (GUT) models one has E/N = 8/3. 
This implies (E/N - 1.93) w 0.75. 

However, in the KSVZ model one finds E/N = 6Q em , where Q em is the 
electric charge of the heavy quark Q. Moreover, in more general KSVZ-type 
models the number N of heavy quarks and their transformation properties under 
SU(3)c are arbitrary. Therefore, as Kaplan || first pointed out, it is possible to 
construct KSVZ-type models with E/N = 2, leading to a very small C ail . In the 
light of the considerable uncertainties of the light-quark ratios w and z, it is even 
possible that C ail = 0. We will see in Chapter [|that KSVZ-type axions with an 
accidentally suppressed photon coupling and masses between 10 and 20 eV are 
astrophysically allowed, particles which we refer to as "hadronic axions." 

2.3.2 Electrons 

Tree-Level Coupling to Electrons: 

In many axion models ordinary electrons carry PQ charges and have fundamental 
Yukawa couplings to the Higgs field. Then the tree-level coupling between axions 



and electrons has the structure of expression ( 2.26|) . The associated coupling 
constant is 

9lT = ^ = C e 0.85 x 10- 10 m eV , (2.42) 

fa 
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Figure 2.3: Radiatively induced axion-electron coupling. 



where C e = X' e /N is a numerical coefficient of order unity. Here, X' e is the shifted 
PQ charge of the electron |lf|, and N the color anomaly. Of course, in KSVZ-type 
models one has C e = 0. 



Radiatively Induced Coupling to Electrons: 



In addition to the tree-level coupling, there exists a loop correction due to the 



axion's interaction with photons (Fig. |2.3|) . Such a radiatively induced coupling 
exists even if X' e = 0. This loop correction can be described by the effective 
coupling constant ffi 



loop 



3a m e 



N \m P 



2(4 + z + w) ln /Aqcd 
3(1 + z + w) \ m e 



(2.43) 



with the cut-off scales f a and Aqcd ~ 150-400 MeV, respectively. 

Considering both the tree-level and the radiatively induced coupling, we ob- 
tain for the total axion-electron interaction Lagrangian 



(2.44) 



2(4 + z + w) 
3(1 + z + w) 



ln 



A 



QCD 



Although the loop correction is of 0(a 2 ) smaller than the tree-level result, it is 
particularly important for KSVZ-type axions because they have g l ™ e = 0. In this 



case, one is left with the radiative coupling (|2.43|) or, numerically, 



9. 



i°° p w 3.4 x 10" 16 



m eV 



^(23.2-lnm eV ) 
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(2.45) 



For hadronic axions with E/N w 2 this is g l °° p ~ 10 14 m e \. 
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2.3.3 Nucleons 



Naturally, if ordinary quarks carry PQ charges they have fundamental Yukawa 
couplings to axions. But in astrophysical considerations one is confronted with 
energies below the QCD scale Aq C d ~ 150-400 MeV. Thus, free quarks do not 
exist and one is interested in the effective coupling of axions to nucleons. This 
interaction arises from two roughly equal contributions. First, there is the tree- 
level coupling of axions to up and down quarks. Moreover, the familiar axion-pion 
mixing provides a second term for the axion-nucleon interaction. The resulting 
effective axion-nucleon interaction can be described by a Lagrangian similar to 
that of Q2.26Q . The axion-nucleon coupling constant is defined as 



9aN = 7 , ( 2 - 46 ) 

J a 

where Cn is the effective PQ charge of the nucleon. This charge is a compound 
parameter combining both parts of the axion-nucleon interaction, i.e. the tree- 
level part and the part due to the axion's mixing with pions. For neutrons and 
protons it is 

C p = (C u -7))Au+(C d -7)z)Ad+(C s -r]w)As , 

C n = (C u - V )Ad+(C d -r}z)Au+(C d -r]w)As , (2.47) 

where r] = (1 + z + w)~ 1 = 0.632 with z = m u jm d and w = m u /m s . Furthermore, 
Aq describes the fraction of the nucleon's spin carried by the quark flavour q. 



Numerically, one finds |16 



An = +0.80 ±0.04 ±0.04, 
Ad = -0.46 ±0.04 ±0.04, 

As = -0.12 ±0.04 ±0.04, (2.48) 

where the first error is statistical, and the second is of systematic nature. 

Let us now consider KSVZ-type axions in more detail. In spite of their van- 
ishing couplings to ordinary quarks, their interaction with nucleons is not zero. 
Actually, the axion-nucleon coupling of KSVZ-type axions is comparable to that 
which appears in models with existing tree-level coupling to quarks, e.g. the 
DFSZ model. The fundamental reason is the generic axion-pion mixing which 
exists even though the tree-level coupling to quarks vanishes. With vanishing PQ 
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charges C u , Cd, and C s of the u-, d-, and s-quarks, respectively, and with ( |2.47| ) 
one obtains explicitly 

C p = -0.34 and C n = 0.01. (2.49) 

These values imply 

9ap = ^ = -5.32xHT 8 m eV , 

J a 

9an = ^^ = 1.57x 10' 9 m eV (2.50) 

J a 

for the KSVZ Yukawa couplings. 
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Chapter 3 

Hadronic Axion Bounds 



Present bounds on the hadronic axions are reviewed; for more general models see 
|J. The possibility of novel hadronic axion bounds, based on neutron star cooling 
and accretion discs around black holes, is deferred to Chapter |6|. Limits on f a (or 
m a ) are usually based on limits on the coupling strength to photons, electrons, 
and nucleons. These limits, in turn, are derived from astrophysical objects like 
the sun or white dwarfs where axion production would increase stellar energy 
losses and thus modify the observed properties of stars ("energy-loss argument"). 
For hadronic axions the relevant astrophysical objects are mainly globular cluster 
stars and SN 1987A. 



3.1 Globular Clusters 

All stars in a globular cluster have nearly the same age and equal chemical proper- 
ties; they differ primarily in one single parameter, their initial mass. Particularly 
interesting for axions are two different kinds of stars, namely those occupying the 
red-giant branch (RGB) and the horizontal branch (HB). The main characteristic 
of the former is their degenerate helium core, whereas the latter are characterized 
by a helium-burning core. In both kinds of stars the core is surrounded by a 
hydrogen burning shell. In the interior of HB and RGB stars axions could be 
produced by different mechanisms, leading to bounds on the axion-photon and 
the axion-electron coupling. 
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3.1.1 Photon Coupling 

In globular-cluster stars photons can be transformed into axions by means of the 
axion-two photon interaction. This effect is referred to as Primakoff conversion. 
It takes place in both the HB and RGB star cores, but, as a consequence of 
the different core densities, Primakoff conversion is much more effective in HB 
stars. Therefore, an additional axionic energy loss would shorten the lifetime 
of HB stars, leading to a reduced HB/RGB number ratio, in contradiction to 



observations. These considerations lead to the bound |17 

£a 77 < 0.6 x 1(T 10 GeV" 1 . (3.1) 

With (pT40D one obtains 

C a77 m e v < 0.3 , (3.2) 

where again m eV = m a /eV and C ail = (E/N — 1.93). Obviously, the limits on 
C Qryy and m a are not independent. In particular, for hadronic axions we have 
C a77 <C 1, and no meaningful limit on m a may be derived. 



3.1.2 Electron Coupling 

Another constraint may be obtained by considering the axionic energy loss via 
the bremsstrahlung process e~ + (A, Z) — > e~ + (A, Z) + a. The consequence of 
this energy loss would be a delay of the helium flash in the red-giant core, resulting 
in a greater core mass at helium ignition. A comparison between the predicted 
and observed core masses at the helium flash yields [fL7| 

g ae <2.5x 10~ 13 . (3.3) 

In the case of hadronic axions, g ae is given by the radiatively induced coupling, 
implying the condition 

m eV (C a77 + 1.3)<35. (3.4) 

As a result of the SN 1987A bound, which will be discussed below, the mass of 
hadronic axions is restricted to the range between 10 and 20 eV. This requires 
C a77 to be less than 0.45-2.2, which can be easily provided in KSVZ-type models. 



3.2 SN 1987A 

Axions could be produced via nucleon-nucleon axion bremsstrahlung N + — ► 
N+N+a in the core of a supernova. As a consequence, the most important bound 
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on hadronic axion masses can be derived from the SN 1987A by considering the 
strength and duration of the observed neutrino signal at the Kamiokande II (KII) 
and the Irvine-Michigan-Brookhaven (1MB) water- Cherenkov detectors. Depend- 
ing on their masses m a , axions either escape freely or are radiated from an axion 
sphere. This energy loss can have three different observable effects: first, it can 
significantly shorten the observed neutrino burst, second, it can cause additional 
counts in the neutrino detectors, and finally, emitted axions could alter the ex- 
tragalactic background light via their radiative decay. 



3.2.1 Nucleon Coupling 

Axions produced in the SN core compete with the standard neutrino cooling 
channel, i.e. they remove energy from the neutrino signal, implying a shortening 
of the observed neutrino signal. The axion luminosity L a depends on the axion- 
nucleon coupling g a ^ and therefore indirectly on the axion mass m a . Assuming 
that the mean free path of axions is larger than the size of the core, axions are able 
to escape freely. In this instance the axion luminosity L a increases with increasing 
axion mass: the greater m a , the stronger g a ?j, and thus the axion emission rate. 
On the other hand, if axions interact too strongly, they are trapped inside the 
core and are emitted from an axion sphere, similar to the well-known concept 



of a neutrino sphere. A simple analytical model [18|, [19] showed that the axion 



luminosity in the trapped regime varies as L a oc m a 16 ^ U , i.e. L a decreases with 
increasing axion mass m a . Naturally, beyond some large coupling g a ^, axions will 
be trapped so effectively that their impact on the SN cooling is again negligible. 
Summarizing these arguments, we can conclude that axions must either couple 
sufficiently weakly or sufficiently strongly to nucleons to avoid a conflict with the 
observed neutrino signals at KII and 1MB. Accurate calculations |2(J |16|, |17]] yield 
the excluded range 

3 x 1(T 10 < g aN < 3 X 1(T 7 . (3.5) 

Assuming a proton fraction of 0.3 inside the SN core, one can calculate an effective 
axion-nucleon coupling C a ^ = y/0.3C% p + 0.7C^ n ~ 0.2, where we have used 
( p.49|) . With g aN = C^m N / f a one finally obtains 



0.01 eV < m a < 10 eV , 
6 x 10 8 MeV < f a < 6 x 10 11 MeV (3.6) 



for the SN 1987A exclusion parameters of hadronic axions. 
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But even axions which are allowed by the energy-loss argument can have mea- 
surable effects, as they might produce additional counts in the water-Cherenkov 
detectors. Engel et al. [P]J showed that such heavy hadronic axions are able to 
induce nuclear excitations in oxygen, a + 16 O — > 16 O*. After that, nuclear deex- 
citations produce 7-rays which can then be detected. Hence, one can estimate a 
range between 

20 eV < m a < 20 keV , 
3 x 10 5 MeV < /„ < 3 x 10 8 MeV (3.7) 

for excluded hadronic axion masses. 



3.2.2 Photon Coupling 

If axions had been emitted from SN 1987A, they might have produced a flux 
of 7-rays due to their radiative decay mode a —>■ 27. Hadronic axions with 
m a ~ 10 eV should have been emitted with a fluence of approximately f a ~ 
5.7 x 10 10 rn~y 2 ^ U cm" 2 (T^|. On condition that the axion-photon coupling g a77 is 
strong enough, the resulting 7-fluence / 7 could have been detected by the gamma- 
ray spectrometer (GRS) on the Solar Maximum Mission (SMM) satellite. The 
fact that the GRS did not detect any signal above the instrument background 
sets an upper limit to the expected 7-fluence / 7 |]T7|, ^2| . 



fy^^C^mf^KOA. (3.8) 

The SN 1987A bound derived in the previous section precludes axion masses 
smaller than 10 eV and greater than 20 eV. Therefore, if we use 20 eV as a 
maximum value for the mass we obtain the conservative limit 

C a77 <7, (3.9) 

which is less stringent than 



3.3 Gamma-Ray Background Limits 

Nontrivial constraints come from the axionic contribution to the extragalactic 
background light. As we will see in Chapter |], hadronic axions could have been 
produced in the early universe. These relic axions could have left their mark due 
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their radiative decay mode a — > 27, provided that the axion-photon coupling g a ~ n 
is sufficiently strong. 

Overduin and Wesson considered the diffuse extragalactic background 
light and searched for axion-decay photons. The absence of a signal implies 

C ary < 0.323, 0.05, 0.015 for axion masses m e y = 5.3, 8.6, 13 . (3.10) 

These limits give strong constraints on the axion-photon coupling, but they can 
be fulfilled assuming an accidental cancellation C ail = 0. 

Furthermore, Ressel |24| searched for photon emission lines in the center of 



three clusters of galaxies, i.e. he considered the flux from a particular region of the 
sky rather than the whole sky. He obtained the bounds C ail < 0.12, 0.059, 0.029, 
0.024, 0.012, and 0.008 for m cV = 3.5, 4.0, 4.5, 5.0, 6.0, and 7.5, respectively. 
However, masses less then 10 eV are already ruled out as a consequence of the 
energy loss argument applied to SN 1987A. 



3.4 Big Bang Nucleosynthesis (BBN) 

The bounds of the previous sections have been frequently discussed so that we 
just gave a brief summary. However, the possible impact of hadronic axions on 
BBN has not been investigated in such detail, justifying a closer look. 

It will be shown in the next chapter that hadronic axions with masses m a = 
10-20 eV would have come into thermal equilibrium after the quark-hadron phase 
transition. At temperature Tp ~ 60 MeV they decouple and are thus highly rel- 
ativistic. The fact that hadronic axions might have an energy density p a compa- 
rable to that of a light neutrino species would affect the outcome of BBN which 
takes place at about 0.05 MeV. During that epoch the universe is radiation domi- 
nated, and the expansion rate H depends on the total energy density p according 
to the Friedmann equation, 

If one takes the existence of relativistic axions at T ~ 1 MeV into account, 
the expansion parameter is increased due to the axion's energy density p a . An 
increase of H means that the all-important neutron-proton ratio n/p, which is 
regulated by the weak interactions, freezes out at an earlier time, when n/p was 
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larger. Therefore, taking hadronic axions into account, more 4 He is synthesized 
than in the standard scenario. 

In order to estimate the influence of thermally produced hadronic axions on 
BBN, one compares their energy density p a with that of a light neutrino species, 
e.g. the electron neutrino. Both particles are relativistic so that 

AJV =feH(i) 4 ' (3 - 12) 

where the factor 4/7 is due to different phase-space distributions and different 
spins of neutrinos and axions. This equation means that the axion's energy 
density is equivalent to an effective number AiV of additional light neutrinos. 
When axions decouple at T D ~ 60 MeV, T a equals T u and thus AiV = 4/7. Below 
Tp the axion temperature scales as T a oc R~ l . The light neutrino decouples later 
at T ~ 1 MeV and g* T 3 R 3 remains constant. In the standard scenario, the total 
number of effective, relativistic spin degrees of freedom does not change until 
BBN takes place, so that T n u oc Therefore, we obtain for AiV at the time 

of BBN 



AAWn = ( — 

Pv 



T - T NS 



4 

-. (3.13) 



This result has now to be compared with the standard model of BBN, something 
that turns out to be anything but simple: In accordance with the standard BBN 
scenario, the mass fraction of helium, conventionally referred to as Y p ( 4 He), de- 
pends not only on A^bbn but also on the uncertainties of the baryon-to-photon 
ratio r\ = 72#/n 7 . The predicted helium abundance must then be compared with 
the observed 4 He abundance, a quantity that is difficult to estimate, implying 
that an upper constraint on AATbbn is afflicted with significant uncertainties. 
Different authors obtain different results: Olive et al. |25[ found the restrictive 
limit of AA^bbn — 0.3. However, Kernan and Sarkar [26| provide the conservative 



bound AA^bbn < 1.53. If we use that value, no danger is ahead as far as hadronic 
axions are concerned because our result ( |3.13| ) is within the above constraint. 



A review of BBN, AA^BBN-limits and related issues can be found in p7| . In a 
nutshell, due to the lack of reliable data, it is premature to infer from BBN-based 
arguments that a hadronic axion with f a ~ 10 6 GeV can not exist. 
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3.5 The Hadronic Axion Window 



In summary, there is a gap of allowed axion masses between 10 and 20 eV. This is 
a result of the fact that the two different SN 1987A exclusion regions — on the one 
hand the "too much energy loss" region ( |3.6| ) , on the other the "too many events in 
detectors" region ( |3.7|) — do not overlap. This 10-20 eV window is open because 
the globular-cluster bounds ( ID] ), (|3.4j ), and the 7-ray background observation 
limits ( |3.10| ) are of no importance if C a77 is accidentally suppressed. According 
to our discussion in Sect. |2.3.1| this is well possible. Consequently, we are not able 
to rule out KSVZ-type axions with a strictly suppressed coupling to photons. 
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Chapter 4 
Relic Axions 



Thermal production of hadronic axions in the early universe is discussed. We 
show that these particles are highly relativistic when they freeze out, implying 
that hadronic axions behave like HDM. Moreover, we estimate the present axion 
density and find that it is comparable with that required in mixed dark matter 
scenarios. 



4.1 Thermal Production in the Early Universe 

If axions were produced in the early universe, they could constitute dark matter. 
In principle, there are two different production mechanisms. First, axions can 
arise as a result of the usual "freeze out process." However, this thermal produc- 
tion is only meaningful if the axion couplings are sufficiently strong, implying that 
m a must be greater than several 10 -2 eV [2S|. Furthermore, there are two non- 
thermal production mechanisms, according to different cosmological scenarios: If 
inflation occurred after the PQ symmetry breaking or if T re h ea t < fa, axions are 
produced by means of the "misalignment mechanism." Before the QCD phase 
transition, the parameter is not at its CP-conserving minimum = 0, but 
somewhere between and ir. Later, at a temperature around T ~ Aqcd, the 
0-field rolls toward = 0, resulting in coherent oscillations which correspond 
to a condensate of zero-momentum axions. On the other hand, if the universe 
did not inflate at all or if inflation occurred before the PQ symmetry breaking 
with T reheat > f a , a network of cosmic axion strings forms at T ~ f a by the Kib- 
ble mechanism, gradually decaying into (massless) axions. At lower temperature 
T ~ Aqcd, these axions acquire a small mass and become nonrelativistic. 
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As a consequence of the axion's Nambu-Goldstone nature, both the misalign- 
ment and the string-decay picture predict that the axion energy density p a is pro- 
portional to m~ 1175 , implying a lower limit on the axion mass around 10 -3 eV. 
Moreover, in the mass range between approximately 10 -5 and 10 -3 eV, nonther- 
mally produced axions would be an ideal candidate for CDM. However, due to 
unknown initial conditions on 0, and uncertainties in the quantitative treatment 
of the string mechanism, the axion mass density p a is not straightforward to 
calculate so that the function Omega(m a ) remains uncertain. 

In this work we are mainly concerned with hadronic axions. For these parti- 
cles, the nonthermal production takes place at a temperature where the axions 
are still in contact with the thermal bath of the universe, implying that nonther- 
mal mechanisms do not take place. Therefore, we will now focus on the thermal 
production of hadronic axions. 

The interactions of axions with the thermal bath are all of the general form 

a + X 1 ++X 2 + X 3 , (4.1) 

where Xi, i =1, 2, 3 are particles of the primordial heat bath. Then, axion 
production is described by the covariant Boltzmann equation, and we may write 
for the number of axions in a co moving volume Y = n a /s 

£ = -^< y - y »>>- (4 - 2) 

where Yeq = n^/s is the equilibrium number of axions in a comoving volume, H 
the expansion rate of the universe, and r a b s = n\{av) is the thermal averaged ax- 
ion absorption rate of the process flj.ip . Furthermore, s = S/R 3 = (2n 2 /45)g*T 3 
is the entropy density, where g* is the total number of effective spin degrees of 
freedom of all relativistic bosons and fermions that are in thermal equilibrium at 
the given temperature. Finally, we used the scaling parameter x = va^jT with 
the nucleon mass and the ambient temperature T. 

Axions decouple from the thermal bath when Y does not change anymore, i.e 



dY/dx ~ 0. According to equation (fO), this freeze out happens when the ratio 



Tabs/ H becomes small, r abs < H. Let us now estimate the decoupling temperature 
Tr, for hadronic axions, i.e. the temperature To at which r a b s ~ H. Under the 
assumption T < Aqcd ~ 300 MeV nucleons already exist in the universe and a 



possible realization of ( |4.1| ) is the axion-pion conversion 

a + N <r-> 7T + N. (4.3) 
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With the couplings g a N = m^/ f a , g^N ~ m^ 1 , and a factor m^ 1 for the nonrela- 
tivistic nucleon propagator one can estimate the interaction cross section for this 
process 

M « (^-— t)' = 7?t (4.4) 

and therefore 

r abs = n N (av) « . (4.5) 

Ja it 

In the case of nonrelativistic nucleons with a chemical potential /i C T the 
number density is given by ~ (m^T) 3 / 2 e~ x . Furthermore, in the radiation 
dominated early universe one has 

^ = 1.67^-^, (4.6) 
mpi x z 

where mpi is the Planck mass. For the relevant temperatures, the particles con- 
tributing to g* are 7, e + , e~, u e , u T , u e , P^, i> T , and a, implying g* = 11.75 so 
that we obtain 

~ 0.60 m N m ™ aT^e"* « 2.7 x 10 6 m e 2 v x^e^ , (4.7) 

with m oV = m a /eV. The decoupling temperature To is reached when this ratio 
equals 1. Fig. |4.1| shows To as a function of hadronic axion masses between 10 eV 
and 20 eV, i.e. for the parameter range of the hadronic axion window. We see 
that the decoupling temperature is approximately 



D 



57-61 MeV . (4.8) 



Although this calculation was only a rough estimate, there is no doubt that 
hadronic axions are highly relativistic when they freeze out. Hence, their equi- 
librium number density is n a = ((3)T 3 /tt 2 so that we find 

^ = 45CP) liB0024 
s 27r 4 g* 

for the equilibrium number per comoving volume. To estimate the present abun- 
dance of relic axions, we return to the Boltzmann equation ( |4.2| ) which can be 
solved easily because Y E q is constant. With the boundary condition Y(0) = 0, 
i.e. initially no axions were present, one finds 



Y{x) = Y EQ 1 - exp 



X -p 

abs 



dx' 



(4.10) 
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Figure 4.1: Decoupling temperature To of hadronic axions as a function of their 
mass m a . 

Strictly speaking, we would have to integrate from the very beginning (x = 0) 
to the present (x = x to day = xo) to find the number of axions in the universe. 



However, today axions are decoupled, and the present ratio T^s/H in (|4.10| ) is 
extremely small. Thus we may integrate forward to x = oo rather then to x = xq. 
Moreover, we are interested in the time after the QCD phase transition which 
occurs roughly at £qcd = ^v/Aqcd ~ 5 so that 



Y(x ) = Y EQ 1 - exp 



QCD 



J- abs 

'x 7 H~ 



dx' 



(4.11) 



Together with expression ( |4.7[ ) and hadronic axion masses m a = O(10 eV), it is 
easy to see that the exponential function exp[— ]'■■■] can be neglected. Hence 
we find that the present number of axions per comoving volume is identical with 
its value at freeze out, i.e. 



Y(x ) = Y EQ (x D ) w 0.024. 



(4.12) 



Now we can calculate the present number density n° a of relic hadronic axions. 
With sq = 2970 cm -3 for the present entropy density, we find 



Y(xq) s ~ 71 cm 



-3 



(4.13) 



One can compare this result with the number density n° = 115 cm 3 of relic 
electron neutrinos and obtains 

n°\ 

w 0.6. (4.14) 

n°J 
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Our calculation of the hadronic axion decoupling temperature ( |4.8| ) and its 
present number density, (|4.13|) or ( |4.14j) , is just an estimate because we have only 



considered the axion-pion conversion process a + N ^ n -\- N. A more careful 
approach includes an accurate computation of all relevant processes. First, these 
are the reactions 

n°n «-> an, 7r p «-> ap, n + n «-> ap, ix~p <-> an. (4-15) 

In addition, there are reactions involving only pions rather than nucleons and 
pions, 

7r + 7r~ <-> a7r°, 7r°7T ± <-> got 111 . (4.16) 

These processes were evaluated by Chang and Choi P0fl . They assert that the pure 
pion processes (|4.16|) dominate over the mixed ones ( 4.15|) for temperatures up to 



150 MeV, which is above our estimated decoupling temperature Tp ~ 60 MeV. 
Therefore, the total interaction rate including all relevant processes is greater than 
our approximation, implying that the ratio Y/H equals 1 at lower temperatures, 
i.e. axion decoupling occurs at temperatures below T D ~ 60 MeV. In fact, Chang 
and Choi found ~ 30-50 MeV. However, axions are relativistic particles 
when they freeze out, and thus the present number density n° of relic axions is 
quite insensitive to the exact value of the decoupling temperature. The reason 
is that n° a in (|4.13|) depends only slightly upon T D , through g*\ x=x ■ However, 



both for Td ~ 60 MeV and Tp ~ 30-50 MeV, the effective number of spin 
degrees of freedom g* is given by 11.75. Therefore, our estimate Q4.12j ) agrees 



with the exact calculation taking all interaction processes into account, which 
yields Y(x ) « 0.02 §. 

4.2 Hadronic Axions as Hot Dark Matter 

The existence of thermally produced hadronic axions, which are highly relativis- 
tic when they freeze out, means that these particles contribute to cosmic hot 
dark matter (HDM). The present hadronic axion number density ( [4.13[ ) can be 



translated into the axion contribution to the present mass density, one obtains 

tt RDM h 2 = 0.007 m eV , (4.17) 

where h is the Hubble constant Hq in units of 100 km sec -1 Mpc -1 . Numerically, 
one finds with h ~ 0.6-0.8 and m a Pd 10-20 eV 

Ohdm « 0.1-0.4. (4.18) 
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This is the right ballpark for a HDM component required in mixed dark matter 
scenarios, pTL |32 
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Chapter 5 

Axions in a Nuclear Medium 



It is the nature of hadronic axions that they couple primarily to nuclear mat- 
ter. Therefore, suitable astrophysical objects for investigating their properties 
are those which consist mainly of nucleons, e.g. supernovae, neutron stars, and 
accretion discs. The impact of axions on SN 1987A was already discussed in 
Chapter while consequences of hadronic axions for neutron star and accretion 
disc physics will be explored in Chapter |6|. The extent to which these objects 
are affected by axions depends strongly on the axion's emission and absorption 
properties. Therefore, we devote the present chapter to the axion emission rates 
and mean free path in neutron stars and accretion discs. 



5.1 Axion Emission in a Nuclear Medium 

5.1.1 Nucleon-Nucleon Axion Bremsstrahlung 

As a consequence of the axion-nucleon coupling, axions can be produced in a 
nuclear medium via the nucleon-nucleon axion bremsstrahlung process N + N — ► 
N + N + a, where N can be a proton p or a neutron n. For this production 
mechanism to be possible, one has to employ a spin-dependent nucleon-nucleon 
potential. An appropriate ansatz is the one-pion exchange (OPE) potential, which 
allows a simple calculation of the bremsstrahlung process. It should be stressed 
that using an OPE potential is just an approximation to the real nucleon inter- 
actions in hot and dense matter. However, since our goal is a simple estimate of 
the axion's influence on accretion discs and neutron stars, the OPE potential is 
a suitable assumption. Furthermore, we assume nonrelativistic nucleons. This is 
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Figure 5.1: One possible Feynman diagram and its exchange graph for nucleon- 
nucleon axion bremsstrahlung. Altogether there are eight graphs, four direct and 
four exchange diagrams. 



justified because the temperatures in neutron stars and accretion discs are around 
100 keV and 10 MeV, respectively. 

In the OPE approximation one has in total eight Feynman diagrams for the 
bremsstrahlung process: four direct diagrams with the axion attached to each 
nucleon line, and four exchange graphs each with N$ 
"pure" processes nn — > nn + a and pp 



iV 4 (Fig. fy). For the 
pp + a with only neutrons and protons, 



respectively, the spin-summed squared matrix element is found to be [33 



3m 2 N 



9aN 



k 2 



mi 



+ 



l 2 



l 2 + m l T[ 

2 



+ 



k 2 l 2 - 3(k - 1) 

(k 2 + m 2 )(l 2 + m 2 



(5.1) 



where g a N = CV tt^n/ fa, N = p, n is the axion-nucleon coupling constant and 
a n = (/2mjv/m vr ) 2 /47r the pion-nucleon "fine-structure constant" with the pion- 
nucleon coupling / ~ 1.05. Furthermore, we have k = p 2 — p 4 and 1 = p 2 — P3 
with the nucleon's momenta p,. For the "mixed" process np — > np + a one obtains 



2 

np 



256vr 2 a 2 (g an + g ap )' 



3m 2 N 



l 2 



4(k-l) ; 



2567r 2 a 2 gl n + g, 



+ 



ap 



3m 2 N 



l 2 + m 2 _ 

1,2 \ 2 



+ 2 



(k 2 + m 2 )(l 2 
l 2 ^ 2 



mi 



k 2 + m 2 / V l 2 + ml 

7T / \ 7T 

k 2 l 2 ~ (k • l) 2 

! (k 2 + m 2 )(l 2 + m 2 )J 



(5.2) 



The larger coefficients compared to (|5 . 1|) are a consequence of the stronger cou- 
plings of charged pions to nucleons. To obtain the axion energy loss rate per unit 
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volume, one has to perform the phase-space integration 




E a Sj2\M\ 2 hf2(l-f 3 )(l-U) 



s 



x (27r) 4 ^(Pi + P 2 - P 3 - Pa ~ P a ), 



(5.3) 



where Pj, j = 1, . . . , 4 are the four- momenta of the external nucleons, P a the one 
of the axion, and fj, j = 1, ... ,4 the phase-space distributions of the nucleons. 
The factor 5 takes identical particles in the initial and final states into account: 
5 = 1/4 for the pure, and 5 = 1 for the mixed processes. The Bose stimulation 
factor for the final-state axion (1 + f a ) was neglected because we will see that 
these particles always escape freely. 

To be able to perform the phase-space integration one needs to specify the 
phase-space distributions fa of the nucleons. For accretion discs and neutron stars 
this is a difficult task. Therefore, as a first approximation we describe the nuclear 
medium in these objects as an ideal gas of protons, neutrons, and electrons, where 
neutrons and protons have approximately the same mass of ph 940 MeV. 
Then, the phase-space occupation function of protons (i = p) and neutrons (i = 
n) is given by the Fermi-Dirac distribution 



As the nucleons are nonrelativistic, it is appropriate to use the nonrelativistic 
kinetic energies Ef m = Ei(pi) — rrii ph pf/2rrii and the nonrelativistic chemical 
potentials //j = /ij — 



erate conditions, limits which pertain to neutron stars and accretion discs, re- 
spectively. For the latter, characteristic temperatures and densities are around 
T < 10 MeV and ps ~ 10 12 g cm -3 , respectively, so that we find for the nucleon 
(nonrelativistic) Fermi energy 

£f,n = P 2 f,n/2 ™n = (3TT 2 Y N n B ) 2/3 /2 rn^lx Y 2/3 MeV < T, (5.5) 

where tib = Nb/V is the baryon number density and = N^/Nb with 
N = n or p the nucleon number fraction. We see that neutrons and protons 
are nondegenerate. In the case of neutron stars, we have T = 100 keV and 
Pb ~ 2p nuc w 5.6 x 10 14 g cm" 3 , where p nuc = 2.8 x xlO 14 g cm -3 is the nuclear 




1 



(5.4) 




The integration in ( |5.3| ) simplifies considerably for degenerate or nondegen 
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density. In Appendix [B| we will show that the proton number fraction Y p is around 

0. 01 so that 

e F , n « 94 MeV > T, e F>p « 4 MeV > T, (5.6) 

1. e. both neutrons and protons are degenerate. 



5.1.2 Nondegenerate Limit 



For a nondegenerate nucleon gas the phase space integration ( |5.3|) can not be 
done analytically. To simplify the squared matrix elements ( |5.1|) and (|5.2j) , we 
first neglect the pion mass in the pion propagator. One then obtains 



with 



~2 , 

9nn — \ 



gin (3 - P) 

9l P (3 - P) 
9an 9ap 



2 

NN 



2567r 2 al _ 



4/9/3) + 



3m^ 



9nn 



2 i 2 



(5 - 2/3/3) 



nn 
pp - 

np 



(5.7) 



> nn + a 
pp + a 

■ np + a . 



Here, g^^ = C]ym N /f a is an effective coupling and (3 = 3((k-l) 2 ). In the 
nondegenerate limit it turns out that (3 = 1.3078 |34l| . Numerically, one finds 
with the axion-nucleon couplings of expression (|2.50|) 



\C 



AT 




nn — > nn + a 
pp pp + a 
np —>■ np + a . 



(5. 



However, neglecting the pion mass will not be appropriate if one considers ac- 
cretion discs: in a nondegenerate thermal medium of nonrelativistic nucleons with 
temperature T, the momenta of the nucleons are approximately given through 

~ 3m nT. Since the temperatures in accretion discs are T w 10 MeV or be- 
low, typical nucleon momenta are |p$| ~ 170 MeV, comparable to the pion mass 
= 135 MeV. Hence, in accretion discs one expects the pion masses to af- 
fect the axion emission rates so that we introduce the "fudge factor" £(T) which 
includes all pion- mass effects. A rough estimate of £(T) in Appendix |A| yields 
approximately 0.5. With a constant matrix element it is then straightforward to 
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perform the integration in expression ( |5.3| ). One obtains for the total energy loss 
rate per unit volume 

^ D = W I ^|f y "^ |A1| «» +y ^ |A1| - +4y » y '"Sl^l». 



it 1 



3? n 2 T' 'I 2 n 2 

—£(T)-?- -^(Y 2 a 2 + Y 2 a 2 + 4Y Y a 2 ) 

= W5 ?(r) ^W 9ND ' <5 - 9) 



where #nd is the total effective axion-nucleon coupling constant for the nonde- 
generate limit. With 1^ »s 0.1, Y n m 0.9, and (|5.8| ) we find numerically 



^ ND = 4.71 x 10~ 8 m a or C^ D = 0.30 . (5.10) 
Finally, we arrive at an axion energy loss rate per unit volume of 

Qa D = 4-89 x 10 27 ergcm-V^v^mfv (5.11) 

and at 

L™ = 9.7 x 10 48 erg b" 1 ^^ m 2 eY . (5.12) 

for the total axion luminosity. M© is the solar mass and M the one of the 
accretion torus. 



5.1.3 Degenerate Limit 

In the degenerate limit, the parameter f3 is zero so that we obtain from ( |5.7| ) the 
squared matrix elements 

J2\^\nn= 2 -^^~9 2 nn (5.13) 



with 



m N 



g 2 an nn^nn + a 

9nn = { 9l P pp^pp + a 

9 2 an + 9 2 ap + gLgL/ 3 np-^np + a. 



lap ' aap&ap 

Numerically, the effective coupling is 





[ 0.01 


nn - 


-> nn + a 


\Cn\ = | 


. 0.34 


pp - 


■> pp + a 




[ 0.338 


np - 


-> np + a 



(5.14) 
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With these matrix elements one can determine the axion volume emission rate. 
In the degenerate limit, there is no need for introducing a fudge factor £(T) 
because the integration in ( |5.3j ) can be performed analytically without neglecting 
the pion masses. It turns out that the = rates must be supplemented with 



a factor [17 



b (u) — 1 arctan 

v ; 6 



u 



u 



3(u 2 + A) 6V2m 2 + 4 
( 2^2u 2 + 4 



x arctan 



it- 



(5.15) 



where m = m n /p F ^ N . For a neutron star with p B « 2p nuc , one has u ~ 0.32 Y N 
ite I 

3n B 



so that we may write F(Y^). Finally, we find for the total emission rate |T7 

1/3 



31 
967680 



7T 



\pp 



I rip 



SItt 5 / 3 ^) 1 / 3 ^ 

3780 m 2 N 
31ir 5 / 3 (3n B )V 3 alT 6 2 

3780 m 2 , 



T 6 Fj/ 3 F(F n ) £ \M\ 2 nn + Y p ^ F(Y P ) £ \M\; 

\ s s 

+ 4Y^F(Y np )J2\M\; 

s 

F(Y n )g an -\-Yp F(Y p )g av + Y n p F(Y np )g np j 

(5.16) 



where go is the total effective axion-nucleon coupling constant for the degenerate 
limit and Y np the effective nucleon fraction for the mixed processes. According 
to Brinkmann and Turner [33], the latter is given by 



yl/3 = _J_ 
np 2V2 



1/2 



2/3 1 



v 2 / 3 V 

J 77, 1 r 



v 2/3 v 2/3 
in T Jp 



(5.17) 



With Y p = 0.01 and Y n = 0.99 one has Y np = 0.06, and we can determine the 
suppression due to non- vanishing pion masses for the nn-, pp-, and np-processes, 
F(Y n ) « 0.64, F(Y P ) « 0.12, and F(F np ) « 0.31. Altogether, we arrive at the 
following effective coupling constant 



# D = 2.04 x 10~° m a or 



0.13. 



(5.1* 



This together with ( |5.16| ) leads to the following expressions for the axion emission 
rate per unit volume, 



QD = 3.74 x 10 29 erg cm" 3 sec" 1 T[ 



MeV m eV 



Pns 



1/3 



(5.19) 
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where m e v = m a /eV and TueV = T/MeV. Finally, we calculate the total axion 
luminosity and obtain 

L^3.4xl0«er gyr - m2v (-|)(^- 2 % 96 , (5 . 20) 
where M is the mass of the neutron star and Tg = T/10 9 K. 



5.2 Axion Absorption in a Nuclear Medium 

The energy loss of a star due to particle emission depends on whether these 
particles are trapped or not. In this section we will show that axions stream 
away freely both from neutron stars and accretion discs. 

The dominant absorption process for axions is inverse axion bremsstrahlung 
N + N + a^N + N. The associated mean free path A a of the axion is related 
to its absorption rate r abs through A" 1 = r abs . One obtains 

A„ — 1 abs — e 



X 



3 — 1 s 

(2tt) 4 5^(P 1 + P 2 -P 3 -P 4 - P a ), (5.2i; 



where E a is the axion energy. The factor e Ea ^ T accounts for the detailed-balance 
relationship which reveals that r em = exp(—E a /T) r a b s , where r em is the axion 
emission rate. A comparison between ( |5.tj| ) and ( |5.21| ) shows that the volume 
emission rate Q a is related to the mean free path by 



i r°° 

Qa = j2^J ^e-^Eadpl (5.22) 



In the case of the nondegenerate limit, one obtains for the axion's inverse 
mean free path [17] 

° " 6mfT'/> i 9 ™' <5 ' 23) 

where x = E a /T. For freely streaming axions, the average energy is determined 
by the nondegenerate nucleon-nucleon bremsstrahlung axion spectrum. We find 
(x) = (E a )/T = 16/7 so that 

A^ D_1 = 3.09 x 10- 9 cm- 1 T M % 5 pj 2 m 2 eY (5.24) 
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or 



Af 3 « 3240 km T° e 5 v P^ 2 ■ (5-25) 



In accretion discs, the maximum temperatures and densities are 10 MeV and 
10 12 g cm -3 , respectively. Then, one obtains for the mean free path of an 20 eV 
axion A a ~ 26 km. The maximum geometrical dimensions of accretion discs are 
about 50 km so that it is not clear if axions escape freely. However, as we have 
used extreme vales for T and p, the average mean free path will be larger than 
26 km so that freely streaming axions are a good approximation. 
In the degenerate case, the inverse mean free path is given by |L7 



iD -i _ ( 3n B )V 3 cer 2 (x 2 + 4tt 2 ) 

247T 7 /3 m ^ (i _ e -x) 



AT = ' V „ ™ ' & (5-26) 



or numerically, 



/ \ -1/3 

A° « 8350 km m e " v 2 T 9 ~ 2 I M. J , (5.27) 

where we have used the average axion energy (x) = (E a ) /T w 3.16, which can be 
calculated using the degenerate nucleon-nucleon bremsstrahlung axion spectrum 
1 17]. For a neutron star with T w 10 9 K and p ~ 2p nuc we obtain 

A° w 6630 km . (5.28) 

Axions of mass m a are not trapped inside a neutron star of radius R if A a > i?. 
For i? 10 km this is fulfilled for 

m a <26eV. (5.29) 

Therefore, axions in the hadronic axion window are allowed to stream freely from 
neutron star interiors. As we will discuss more precisely in Chapter ^, one expects 
parts of the neutron-star matter to be in a superfluid state. If so, the mean free 
path (|5.28|) is significantly enhanced. 
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Chapter 6 

Impact on Black Hole Accretion 
Discs and Neutron Stars 

We explore how the emission of axions with parameters in the hadronic axion 
window affects the evolution of gamma-ray bursts. We assume that a subclass 
of short bursts with a duration of 0.1-1 seconds is driven by the accretion of hot 
plasma by a black hole. Axion production in this accretion torus would provide 
a significant energy loss, altering the final 7-fluence. Moreover, we investigate 
the impact of an hadronic-axion energy loss on the cooling timescale of neutron 
stars. We find a relation between their present surface temperature and the axion 
mass. Our results are then compared with the observational data of the pulsars 
PSR 1055-52, PSR 0630+178 (Gemina), and PSR 0656+14. Finally, we discuss 
effects of nucleon superfluidity on our calculation. 

6.1 Gamma- Ray Bursts 
6.1.1 General Picture 

Gamma-ray bursts are short and intense bursts of photons with energies in the 
range between approximately 100 keV and 1 MeV. In fact, gamma-ray bursts are 
the electromagnetically most luminous objects observed in the universe. Since 
their discovery in the late sixties, several satellites have observed them, and nu- 
merous theories have been developed to explain their physical nature. A milestone 
in the exploration of gamma-ray bursts was a remarkable finding of the BATSE 
detector, which was launched in the spring of of 1991. BATSE measured a per- 
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fectly isotropic angular distribution of gamma-ray bursts and a deficiency of faint 
bursts p5| . This led to the suggestion that gamma-ray bursts are cosmological, 
an assumption that was confirmed by the discovery of x-ray transient counter- 
parts by Beppo-SAX [|36[| , and later by the discovery of optical |37j and radio 
counterparts fl38f . For the first time, it was possible to rule out a wide class of 
gamma-ray burst theories! 

In spite of these discoveries, the origin of gamma-ray bursts remains unclear. 
However, there exists a promising generic model which is in agreement with obser- 
vation: the "fireball model," where an ultra-relativistic energy flow is converted 
into radiation. The initial energy flow can be an electromagnetic Pointing flux 
or kinetic energy of highly-relativistic particles. As the conversion in radiation 
takes place in an optically thin region, gamma-ray burst spectra are nonthermal, 
in agreement with observation. 



6.1.2 The BHAD Model 

The nature of the central engine that drives the gamma-ray burst remains unclear 
because the "inner engine" that powers the relativistic fireball is hidden from 
direct observations. While many models have been proposed, those currently 
favored are all based on the rapid accretion of matter by a black hole. These 
models are referred to as black hole accretion disk (BHAD) models. It must be 
stressed that BHAD models are only able to explain short gamma-ray bursts with 
durations in the range 0.1-1 s. Several promising candidates for the progenitor 
system of a BHAD have been proposed, the most popular being neutron star 
binaries, which we presently focus on. 

The merger of two neutron stars leads to a hot and dense central object with 
masses around 3 M , surrounded by a toroidal cloud of hot plasma gas with a 
mass of about 0.1-0.2 M Q . Within a timescale of a few milliseconds, the 3 M Q 
core will most likely collapse to a black hole. However, a significant amount of 
mass, typically 0.03-0.3 M , has too much angular momentum to be swallowed 
immediately by the black hole, and an accretion disc forms. Such accretion discs 
have masses between several 10~ 2 M Q and a few 10 _1 M Q . Their densities are 
around 10 10 -10 12 g cm -3 , their temperatures in the range between 3 and 10 MeV. 

Within this accretion torus neutrinos of all flavors are produced by different 
processes. The bulk of the neutrino luminosity is provided by the electron neu- 
trinos v e and z/ e , which are mainly produced by the charged-current processes of 
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electron and positron capture on protons and neutrons p9|, 



e + p — > n + v e , 

e + + n — > p + v e . (6.1) 

Second is the muon and tau neutrino production, principally via electron-positron 
pair annihilation 

e~ + e + — > Vi + i>i , i — [A or r . (6.2) 

Depending on the specific BHAD model, electron neutrinos are emitted from a 
neutrinosphere with an outer radius of 70-100 km. Moreover, muon and tau 
neutrinos may be trapped, depending on the thermodynamic properties of the 
torus 0. The total neutrino luminosities L v are typically in the ballpark of 
10 53 erg s _1 . 

Now we come to the crucial point concerning gamma-ray bursts. The emitted 
neutrinos and antineutrinos interact with each other and annihilate into electron- 
positron pairs. The efficiency of this energy deposition E uD in form of an electron- 
positron plasma is usually expressed in terms of the parameter q uD = E uP /L u , 
where E v9 is the rate of energy deposition by neutrino-antineutrino annihilation. 
The efficiency q ui> depends on the number densities and energies of the initial 
neutrinos. According to 0, typical values for q vV are around 0.004, i.e. only a 
half percent of the initial neutrino energy is transformed into thermal energy of 
an electron-positron plasma, the fireball. Integrated over the accretion timescale 
t pa 5-50 ms, the total energy deposition is of order 10 49 erg. Moreover, it is 
possible that the fireball does not expand isotropically, but is beamed into axial 
jets with a fraction /q = 25f2/47r of the whole sky. Including a beaming factor 
/n ~ 0.1-0.01, it follows that the BHAD scenario is able to account for short 
gamma-ray bursts with energies in the range of 10 50 -10 51 erg. 



6.1.3 Influence of Hadronic Axion Emission 

Let us now examine to what extent hadronic axions could alter the BHAD model 
for gamma-ray bursts. During the process of accretion, axions could be pro- 
duced by nucleon-nucleon axion bremsstrahlung. Clearly, in order to maintain 
the BHAD model, the total axion luminosity L a may not exceed the neutrino 
luminosity L u , which is of about 10 53 erg s _1 . 

We consider two BHAD models that have been investigated in great detail [0 . 
The first describes the attractive force of the central black hole by a Newtonian 
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Table 6.1: Temperature T, density p, mass M and neutrino luminosity L v of an 
accretion disc for the Newtonian and the Paczyhski-Wiita model. T, p, M, and 
L v are the associated averaged quantities. 





Newtonian model 


Paczyhski-Wiita model 


T [MeV] 


2-10 


2-10 


T 


6 


4 


p [10 12 g cm -3 ] 


0.01-1 


0.01-0.3 


P 


0.3 


0.1 


M[M ] 


0.16-0.29 


0.002-0.031 


M 


0.23 


0.02 


L v [10 52 erg s- 1 ] 


1-12 


0.14-6.7 


Lv 


5.6 


4.4 



gravitational potential, which is known to be singular at the origin r = 0. The 
second takes general-relativistic effects into account by allowing one to reproduce 
the existence of a last stable circular orbit at a radius of 3r s = 6GM^/c 2 , where 
r s is the Schwarzschild radius. This model is referred to as the Paczyhski-Wiita 
model. The resulting temperatures T, densities p, and accretion disc masses 
M for both models are shown in Tab. |TT] as well as the mean values of these 
quantities. One sees that the torus masses in the Paczyhski-Wiita model are 
much smaller than those of the Newtonian model. This is a result of the very 
strong gravitational potential of the Paczyhski-Wiita model, implying that the 
accretion rates, i.e. the mass accretion into the black hole per unit time, of this 
model are much greater than those in the Newtonian description. Therefore, the 
torus mass of the latter is significantly higher than the one of the former model. 
However, it might be that the Paczyhski-Wiita model underestimates the torus 
mass. 

Now, we want to estimate the axion luminosity L a . As we have seen in the 
previous chapter, axions stream away freely from accretion discs, implying that 
axion volume emission takes place. The corresponding energy loss is given by 
expression ( |5.12| ). It is linear in the mass M and the density p so that one may 



use the average values given in Tab. |6.1| . Then, one finds for the two models 

6.7 x 10 47 erg s _1 Tj^y m lv Newtonian model ^ 
1.9 x 10 46 erg s -1 T^ v ml v Paczyhski-Wiita model. 
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Newtonian Model 



Paczynski-Wiita Model 




6789 10 6789 10 

T [MeV] T [MeV] 



Figure 6.1: Axion luminosity relative to the one of neutrinos L a /L u as a function 
of the temperature T and the axion mass m a . The left and right panels are based 
on the Newtonian and Paczynski-Wiita model, respectively. 



Now, these luminosities have to be compared with those of the neutrinos given in 
Tab. 6J] . This is graphically done in Fig. |6.1| , where the ratio L a /L u as a function 



of the temperature T and the axion mass m a is plotted. 

Obviously, in the case of the Newtonian model, hadronic axions could alter 
the existence of gamma-ray bursts, whereas for the Paczynski-Wiita model the 
axion luminosities are at least one order of magnitude smaller than the one of the 
neutrinos. This is a consequence of the small torus masses in the Paczynski-Wiita 
model. 

Finally, we want to estimate the maximum axion production rate which can 
be obtained for the maximum temperature T « 10 MeV. With (|6.3|) we find 



. 2.1 x 10 51 ergs 1 m^y Newtonian model (6 4) 

6.0 x 10 49 erg s _1 mg V Paczynski-Wiita model. 

As a result, one obtains an upper bound on the allowed axion mass of approxi- 
mately 

{5 eV Newtonian model 
(6.5) 
27 eV Paczynski-Wiita model. 

It is not clear which of both models gives an appropriate description of the 
accretion process, or whether they are a good approximation of the black hole 
gravitational potential at all. Therefore, we estimate that axions in the parameter 
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range between 10 and 100 eV are in principle able to alter the evolution of gamma- 
ray bursts significantly. Note that this range includes the hadronic axion window 
of 10-20 eV, i.e. the existence of hadronic axions would be relevant for gamma-ray 
burst physics. 

Presently, nobody knows whether the BHAD model gives the correct descrip- 
tion of gamma-ray bursts. However, if this picture proves to be correct, at least 
for short gamma-ray bursts, one could calculate reliable bounds on the axion 
mass. On the other hand, if hadronic axions were detected by the solar-axion 
search experiments |4(], 4T[] , the BHAD model certainly would have to include 
axion emission in a detailed description. 



6.2 Neutron Stars 

6.2.1 A Simple Model for Axion Cooling 



We showed in the Sect. |5.2| that axions in the hadronic axion window are allowed 
to stream freely from neutron star interiors, implying an additional energy loss 
which could have observable consequences which shall now be explored. There- 
fore, we set up a simple model of late-time neutron star cooling, assuming that the 
star cools via axion volume emission. Furthermore, we neglect all other cooling 
mechanisms because axions only alter neutron star cooling significantly if they 
are the dominant cooling channel. The corresponding axion luminosity was given 
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L » = 3 ' 4Xl049 ^ yr " lm -(^)(^)" 2/ ^- (6-6) 

This result is based on the assumption that the nuclear medium can be treated 
as an ideal, degenerate Fermi gas of neutrons and protons. However, the real 
nuclear equation of state of the superdense matter inside a neutron star is still 
unknown. According to the BCS theory of superfluidity, one expects parts of the 
neutron star to be in a superfluid state. At lower densities p ~ 10 13 -10 14 g cm -3 
particles are expected to be in a 1 So superfluid state when the temperature falls 
below the critical temperature T cr . A second superfluid state 3 P 2 occurs in more 
dense matter with p ph 10 14 -10 15 g cm" 3 , implying that neutrons in the core of 
neutron stars should be in this superfluid state. As the density of protons in a 
neutron star is much smaller than the one of the neutrons, they are expected to 
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be in the 1 So state. However, superfluidity in neutron stars is a difficult matter, 
and many of its aspects are not yet understood. Particularly controversial is the 
proton 1 So pairing. But, if superfluidity occurs, it has different effects on the 
cooling of neutron stars; the most important in the context with axion emission 
is the suppression of our estimated axion- nucleon interaction (p. 18 ). Moreover, 
strong superfluidity reduces the heat capacity C of the nucleon gas. We include 
these effects by introducing an effective axion luminosity and heat capacity, 

Lf=^L a = ^ATl and C cS = -j- C, (6.7) 

Jem Jem JC 

respectively, where f em and fc are fudge factors. 

The effective axion luminosity is linked with an energy loss of the neutron 
star. According to our simple model of nuclear matter, the total, unsuppressed 
heat capacity C of a neutron star is that of a degenerate ideal gas of neutrons, 
protons, and electrons, 

7r 2 / N n N„ N e \ . . 

C = —T I — — H — H I , (6.8) 

where N n , N p , and N e are the numbers of neutrons, protons, and electrons and 
e Fn , e Fp , and ep ne are the related Fermi energies. For nonrelativistic nucleons 
i = n,p one has p F>i = (3^Y inB )^ = p F , N Y^ 3 and e Fti = p%/2m N = e F>N Y?" 
whereas e F>e = p Ffi for the relativistic electrons. Charge neutrality implies n e = 
n p and thus p F>e = p F>p so that e F)£ = p FtP . Finally we obtain 

7T 2 N B T ( 1/3 1/3 2/3 £ F>N 

9c- \ n P P r, 

* e F,N V PF,N 

= —^^(0.997 + 0.215 + 0.009), (6.9) 

2 £f,n 

where we have used Y p = 0.01. Obviously, neutrons, protons, and electrons 
contribute around 82%, 17% and 1%, respectively, to the total heat capacity 
C. If we express C in terms of the solar mass and the nuclear density, we find 
explicitly 

^, 42xl0 ^T„(|-)(^)- 2/3 , B T, (6,0) 

As axions are considered to be the dominant cooling channel, one obtains as a 
result of energy conservation 

dE NS dT ff 

-W- c ^Tt-~ L ^ (6 ' n) 
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where E^s is the total thermal energy of the neutron star. With ( |6.7| ) and ( |6.10| ) 
we obtain the equation 

> 9 § = -liT 9 6 . (6.12) 

JC «£ Jem 

This equation can be easily integrated with the result 

At = t f -t t = ^ T 9 (/r 4 (l - , (6.13) 

where T 9 (i) = Tj/10 9 K and T 9 (f) = Tj/10 9 K are the internal temperatures 
of the neutron star, taken at the initial and final times ti and tf, respectively. 
Furthermore, we have introduced the total fudge factor / sup = f cm / f c . Without 
introducing any large errors, T 9 (f)/T 9 (i) can be taken to be zero. Then ( |6.13|) 
simplifies to 



t yI « 0.010 m e - v 2 /su P T 9 (/)- 4 (6.14) 

or 

v,... = n 39 , , 

' sup eV yr 



T C9 = 0.32 fg m e " v 1/2 t-V 4 , (6.15) 

where t yT = tf/yx is the neutron star age in years and Tc 9 = T 9 (f). It should be 
stressed that T C9 is effectively the core temperature because we do not consider 
axion production in the neutron star crust. Consequently, equation (|6.15f) allows 
us to calculate the internal temperature Tc of a neutron star depending on its 
age t yT , the axion mass m a , and the fudge factor / sup . A remarkable feature of 
the solution ( |6.15| ) is that it depends neither on the mass M nor on the density 



p of the neutron star. 

In order to compare our prediction with observational data, it is appropriate to 
rewrite equation ( 6.15| ) in terms of the observable surface temperature T$, instead 



of the internal temperature Tq. Strictly speaking, the temperature detected by a 
distant observer is not T$, but the red-shifted temperature T|° = T s ^/l — R g /R, 
where R g is the gravitational radius. For neutron stars, the stellar radius R is 
only 2-3 times larger than R g , i.e. R/R g = 2.5. Therefore, we obtain the relation 

Ts^l.29T~. (6.16) 

Furthermore, as a consequence of thermal conductivity, the core temperature Tc 
is related to the surface temperature T s . Approximately, the following equation 
holds |2j: 

T C9 = 0.128 MM , (6.17) 
\9uJ 
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Table 6.2: Measured effective surface temperatures T|^ bs and derived bounds 
on the ratio / SU p/ m eV' wnere /sup is a fudge factor which takes the suppression of 
axion emission due to superfluidity into account. 



Source 


lOgtyr 


logT 5 - obs [K] 


Exclusion range for f SVLp /ml v 


PSR 1055-52 


5.73 


5.84-5.91 


/su P /m c 2 v <1167 


Gemina 


5.53 


5.67-5.80 


/sup/m c 2 v <43 


PSR 0656+14 


5.00 


5.93-5.98 


/sup/m c 2 v <975 



where Tse = T5/IO 6 K and #14 is the gravitational acceleration on the neutron 
star surface in units of 10 14 cm sec -2 . Eliminating the temperature Tc9 between 
the equations (|6.15|) and ( |6.17|) , leads together with (|6.16|) to 



= 1.282 gli 4 f snp ^° #V«> • (6.18) 

The gravitational surface acceleration is typically in the range of gu ~ 2-3. 
Therefore, with g u = 2.5 and the logarithmic scales logTg and logt yr , we finally 
arrive at 

where log is the logarithm to base 10. Now we are in a position to compare our 
predictions with the observational data from ROSAT. Axions with masses m a 
are execluded if the predicted temperatures log T|° of expression (|6.19|) are below 



logT~ « 0.138 log ( ^f- ) - 0.138 logt yr + 6.208, (6.19) 



those observed, i.e. T§° < Tg° ohs . The measured effective surface temperatures 
T™ hs and the ages of three pulsars PSR 1055-52 |43|, PSR 0630+178 (Gemina) 
44j], and PSR 0656+14 j|5|| are listed in Tab. |6.2| . Additionally, in the fourth 



column the exclusion ranges for /sup/^ev are gi yen f° r eac h source. To rule out 
the whole hadronic axion window between 10 and 20 eV, we obtain maximal 
allowed fudge factors / sup of order 10 5 , 4 x 10 3 , and 10 5 for the pulsars PSR 1055- 



52, Gemina, and PSR 0656+14, respectively. Fig. |0| shows the predicted cooling 
curve of a neutron star according to our model. The axion mass is m a = 10 eV 
and m a = 20 eV for the dashed and solid curves, respectively. Furthermore, the 
fudge factor was supposed to be / sup = 10 3 . In spite of this extreme suppression 
of the axion emission, the associated axion cooling predicts present-day surface 
temperatures T£° ohs which are far below the measured data, i.e. the observed 
neutron stars are too hot to have been cooled by hadronic axion emission. 
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Figure 6.2: Red-shifted surface temperature T|^ bs versus age for axion cooling. 
The dashed and solid curves correspond to axion masses of 10 eV and 20 eV, 
respectively. The fudge factor / sup was assumed to be 10 3 . 

6.2.2 Numerical Cooling Calculations Including Axions 

Recently, Umeda et al. || have investigated the thermal evolution of neutron stars 
including axionic energy losses. They intended to compete with the low mass SN 
1987A limit of equation Q3.6| ), which implies the upper limit m a < 0.01 eV on the 
axion mass. Just as we did above, they assumed axion emission due to nucleon- 
nucleon bremsstrahlung. They based their calculation on three equations of state, 
a stiff, a medium, and a soft one; corresponding to different compressibilities of 
the nuclear matter. Furthermore, they took 3 P 2 Cooper pairing of neutrons into 
account. Umeda et al. compared their results with the data of the PSR 0656+14 
pulsar and obtained the limits 

{0.33 eV stiff model 
0.08 eV medium model (6.20) 
0.06 eV soft model, 

which are less stringent than the one of SN 1987A . 

Nevertheless, we can gain from their calculation because it has implications 
for the hadronic axion window. Umeda et al. assumed that axions do not get 
trapped in the interior of neutron stars. Hence, their constraint applies as long 
as the free streaming condition A a > R is satisfied. However, we have seen in 
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Sect. [5.2| that hadronic axions stream away freely, i.e. the result of Umeda et al. 
excludes axions with parameters of the hadronic axion window. 

Over and above that, we can use their limits to estimate of the fudge factor 



/sup. With ( |O0D and the condition f SU p/m 2 eV > 975 for the PSR 0656+14 we find 




stiff model 

medium model (6.21) 
soft model. 

The fudge factor / sup is a compound parameter taking the suppression of axion 
emission and heat capacity due to superfluidity into account. In the case of the 
stiff equation of state, neutrons are assumed to be strongly superfluid, implying 
a significant reduction of their heat capacity — actually it becomes smaller than 
the one of the electrons so that neutrons effectively do not contribute to the 
total heat capacity C of the neutron star. We have seen above that neutrons 
contribute around 80% to the total heat capacity. Therefore, for the stiff model 
we estimate a heat-capacity fudge factor fc of about 5, implying f em « 530. In 
models with a softer equation of state, the reduction of the heat capacity as a 
result of superfluidity is smaller, and the fudge factor is difficult to determine. 
Approximately, we use fc ~ 2 and fc ~ 1 for the medium and the soft model, 
respectively. Altogether we have 

stiff model 

fem ~ { 12 medium model (6.22) 
soft model 

for the suppression of the axion emission as a consequence of neutron superfluidity. 

Now, we are in a position to specify the axion trapping inside a neutron star 
more precisely. The mean free path A^ was given in equation (|5.28| ), however, 
without taking superfluidity into account. The effective mean free path is given 
by Af = / em A°, i.e. 

Af « 6630 km f em m c ~ v 2 . (6.23) 

Axions get trapped if A° ff is comparable to the neutron-star radius R ~ 10 km. 
Therefore, from ( |6.22| ) we obtain the following conditions for free streaming ax- 
ions, 

' 593 eV stiff model 
m a < ^ 89 eV medium model (6.24) 
51 eV soft model. 
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Hence, we conclude that axions do not get trapped in the interior of neutron stars 
if their masses are less than 50-600 eV. 



6.3 Is the Hadronic Axion Window Closed? 

We have re-examined the possibility of an astrophysically allowed KSVZ-type 
axion which has a strongly suppressed coupling to photons, and we have con- 
firmed that such "hadronic axions" with masses between 10 and 20 eV are not 
excluded by previous arguments. Axions in this window were thermally produced 
in the early universe, implying that they constitute HDM with Q a m 0.1-0.4, in 
agreement with the recent work of Moroi and Murayama ||. Hence, hadronic 
axions not only solve the strong CP problem, but they are also of cosmological 
importance. 

As hadronic axions mainly couple to nucleons, we have investigated their im- 
pact on accretion discs and isolated neutron stars. We have found that hadronic- 
axion production in accretion discs would have observable consequences if these 
objects really provide an explanation for short gamma-ray bursts as in the BHAD 
model. This model is based on coalescing compact binaries such as two neutron 
stars, implying the emission of gravitational waves which could in future be de- 



tected by the laser interferometers LIGO, VIRGO, TAMA, and GEO |§. If 
so, a time-coincident observation of gamma-ray bursts and gravitational waves 
would be strong evidence for the BHAD model. However, even if the BHAD 
model proves to be true, it is premature to rule out hadronic axions because 
their emission rate from the accretion torus depends significantly on details of 
the BHAD model. On the other hand, a future detection of solar hadronic axions 



via resonant absorption in Fe [40, 41] would have important consequences for 



the BHAD scenario, perhaps even calling this explanation for gamma-ray bursts 
into question. 

The most significant limits on hadronic axions arise from old neutron stars. 
Our simple cooling model indicates that hadronic axions accelerate the cooling 
process of these objects significantly, with the consequence that neutron stars 
today are actually too hot to be in accord with hadronic axions. Even in the 
presence of superfluidity, which suppresses the emission rates, this discrepancy 
remains significant. The same conclusion is reached on the basis of the numerical 
neutron-star cooling simulations by Umeda et al. ||. Conversely, if the existence 
of hadronic axions were verified, the standard cooling mechanism for neutron stars 
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would have a serious blemish. In this case, internal heating effects would have 



to be included |47]] or the interpretation of the observed soft x-ray components 
of the pulsars PSR 1055-52, Gemina, and PSR 0656+14 |3|, g§ 
blackbody spectra would have to be reconsidered. 



as thermal 
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Appendix A 



Pion Mass Effects in the 
Bremsstrahlung Process 



It was mentioned in Chapter || that the influence of non-zero pion masses m n on 
the axion emission via nucleon-nucleon axion bremsstrahlung is a temperature 
dependent effect, and that the necessary phase-space integrations can not be done 
analytically in the case of the nondegenerate limit. It is therefore appropriate to 
estimate these effects as a function of the temperature. Recall that the squared 
matrix element for nucleon-nucleon axion bremsstrahlung is 



m 



N 



k 2 + m 2 

k 2 1 2 



+ 



l 2 



+ 



l 2 + m 2 
3(k-l) 2 



(k 2 + m 2 )(l 2 + m 2 )J ' 



where k = p 2 — P4 and 1 = p 2 — P3 with the nucleon's momenta p« 



For zero pion masses (|A.l ) becomes 



trii=0 



m 



(A.l) 
1...4. 

(A.2) 



N 



with (3 = 3((k-l) 2 ) = 1.3078. 

In a thermal medium the momenta of nonrelativistic nucleons are on average 
(p 2 ) ~ 3mArT, so that one can estimate k 2 = l 2 « QtunT. Then, we find for the 
momentum dependent terms of ( [A.l ) 

2 

= m (a.3) 



k 2 + ml 



l 2 



\2 + m 2 



1 



1 + 



6m N T . 
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Figure A.l: The factor £(T) as defined in ( |A.3|) . 



and 



k 2 l 2 - 3(k 1) ; 



. «^(T)-/3e(r)- 

(k^ + m£)(F + m~) 
Altogether, we finally obtain for the squared matrix element 

256tt 2 g 2 aN a 2 n 



m 



: £(T) (3 - /3) 



A? 



(A.4) 



(A.5) 



(A.6) 



i.e. all effects attributed to a non vanishing pion mass m w are combined in the 
overall factor £(T). In Fig. [A.l| , the function £(T) is plotted for temperatures 
between and 40 MeV. This simple estimation of pion mass effects is in good 
agreement with the exact calculation of Raffelt and Seckel [34" . However, our 
approximation underestimates the reduction of axion emission rates due to the 
pions mass: At temperatures around 10 MeV the exact result predicts a reduction 
of about 50%, in contrast to our estimation, which yields 57%. 
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Appendix B 

Nucleons and Electrons in 
Neutron Stars 



We consider a nuclear medium which is made up of protons, neutrons, electrons 
and neutrinos. We assume that all particle species are in chemical equilibrium 
with each other as result of the reaction e + p <-> n + u e . In addition, neutrinos 
are allowed to stream freely away, i.e. their chemical potential vanishes //„ = 0. 
The phase-space occupation function of each species i — p,n, e, u e is given by the 
Fermi-Dirac distribution 

exp[( J B i (p i ) - Hi)/T\ + 1 

with the energy Ei(pi) 2 = mf + p 2 and the chemical potential \x%. The value of 
fj,i at T — defines the Fermi energy E\ i = mf + p\ i in terms of the Fermi mo- 
mentum pp,i, which is related to the particle's number density by n« = 
with g — 1 for neutrinos and g — 2 for nucleons and electrons. For a nonrelativis- 
tic particle species (neutrons and protons in neutron stars), it is appropriate to 
use the nonrelativistic Energie Ef m = Ei(pi) — rrii ~ p 2 /2 i — n,p and the 
nonrelativistic chemical potential \Xi = //j — to,. 

It is useful to introduce the neutron and proton number fractions Y n and Y p , 
respectively. They are defined as 

y n = n n /n B and Y p = n p /n B , (B.2) 

where n# = (N p + N n )/V is the total baryon number density. Consequently, 
Y p + Y n — 1. Furthermore, the electron fraction Y e = n e /riB describes the number 
of electrons per baryon. As a result of charge neutrality, the electron fraction is 
equal to the proton fraction, Y e = Y p . 



59 



The number fractions Y p and Y n can be determined easily for neutron stars 
if one works on the assumption that the star is transparent to neutrinos, i.e. 
fi Ue = 0. This is justified since we cosider quite cold neutron stars. One then 



obtains the proton fraction in the case of a degenerate medium [4c 



with 



n p _ 1 fl + 4Q/m n x 2 n + A{Q 2 -m 2 e )/m 2 n x 4 ^ 3/2 



, , , . (B.4) 

n n 8 V 1 + l/xl 1 

where Q = m n —m p = 1.29 MeV, x n = pF,n/m n , andpF, n is the Fermi momentum 
of the neutron. From that, the neutron and electron number fractions follow 
immediately through Y n = 1 — Y p and Y e = Y p . 

We are concerned with neutron stars which have characteristic densities of 
Pb ~ 2p nuc w 5.6 x 10 14 g cm -3 and core temperatures less than 100 keV. Then, 
one explicitly finds Y p w 0.009, i.e. nuclear matter in neutron stars is approxi- 
mately made up of 99% neutrons and 1% protons. However, it should be stressed 
that these values depeds strongly on the underlying equation of state. Moreover, 
a neutron star is not just a homogeneous ball but consists diverent layers like the 
outer crust, inner crust, outer core, and inner core. Another problem is that in 
a dense medium nucleon interactions become important so that the ideal gas is 
no longer a good approximation of the equation of state. As a consequence, the 
nucleon effective mass m* N is less than its vacuum value, and its thermal excita- 
tions depend on a uncertain dispersion relation. Summing up, we conclude that 
the proton fraction in neutron star is a few percent of the neutron number. 
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